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Exotic excitations arise at the interface between a three-dimensional topological insulator (TI) 
and superconductors. For example, Majorana fermions with a linear dispersion, E ~ k, exist in a 
short it Josephson junction on the TI surface. We show that in these systems, the Andreev bound 
states spectrum becomes nearly flat at zero energy when the chemical potential is sufficiently away 
from the Dirac point. The flat dispersion is well approximated by E ~ k , where N scales with 
the chemical potential. Similar evolution from linear to flat dispersion also occurs for the subgap 
spectrum of a periodic superconducting proximity structure, such as a TI surface in contact with a 
stripe superconductor. 
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Moving at "the speed of light", vp, massless Dirac 
electrons on the surface of a three-dimensional Z 2 topo- 
logical insulator (TI) can not be localized by scattering 
from nonmagnetic impurities [TJ [5] , nor can they be eas- 
ily confined by electrostatic potentials due to Klein tun- 
neling Proximity coupling to ferromagnetic or su- 
perconducting order can however open up a gap in the 
spectrum, thus rendering excitations massive [TJ [5] • An 
intriguing possibility is to engineer new massless excita- 
tions by confining and coherently mixing Dirac electrons 
and holes using two or more superconductors with def- 
inite phase difference [4]. For example, Fu and Kane 
showed that a Josephson junction on the surface of a 
TI with a phase bias of it is a one-dimensional quantum 
wire for Majorana fermions, which can be further manip- 
ulated by using tri-junctions [4 . Signatures of Majorana 
fermions in such structures have been reported in recent 
experiments |6] . 

In this Letter, we demonstrate a drastically different 
regime for the same, albeit slightly more general, Joseph- 
son structures considered by Fu and Kane. This regime 
features massless zero energy excitations that are al- 
most dispersionless, i.e. with vanishing group velocity 
(dE/dk = 0).We elucidate the scattering kinematics be- 
hind the nearly flat dispersion at zero energy using simple 
models, and verify the results with self-consistent calcu- 
lations. We find it striking that in such simple structures, 
which are now available in experiments, the low energy 
excitation can be easily tuned all the way from E ~ k 
to E ~ k , where N is large, by increasing the chemi- 
cal potential. By extending such junctions into a class of 
periodic superconductor- TI proximity structures, we fur- 
ther show that these states become a flat band near zero 
energy. 

The Josephson junction is schematically shown in 
Fig. |TJl) . Two s-wave superconductors are patterned 
on the TI surface. Due to the proximity effect, the 
S-TI interface becomes a 2D superconductor (S). The 
S-TI-S junction can be well described by the following 
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FIG. 1: (color online) a) Schematic of a Josephson junction 
on the surface of a topological insulator (TI). The two super- 
conducting leads (S) have a phase difference tt. A is the su- 
perconducting gap, and w is the junction width (not to scale), 
b) Specular Andreev reflection in the regime E > /i. c) Retro- 
reflection for E < ii. d) Dark lines show the (k y , fi) values for 
the zero energy Andreev bound states for w = IOSuf/A and 
L — > ao. k y is in unit of A/Fivf- 

Bogoliubov-Dirac Hamiltonian introduced in Ref. [4], 

T-L = hv F {cr x ky + iT z a y d x ) + T z fi(x) + T y a y A(x). (1) 

Here Tj (<7j) are the Pauli matrices in the particle-hole 
(spin) space. The system is translationally invariant in 
the y direction, and k y is the momentum along y. In 
the TI region of length w, the superconducting order pa- 
rameter A(x)=0, while it is constant A deep into the 
superconductor. The chemical potential \i can be tuned 
by applying a gate voltage. In general, its value can differ 
in the TI and S region, but for simplicity, we assume it 
is uniform in all regions. Also, we will focus on the case 
of phase difference of 7r across the junction. 

We first give a heuristic argument for the existence 
of two regimes. A Dirac electron in the TI region in- 
cident on S will be Andreev reflected into a hole if its 
energy is below the superconducting gap (E < A). In 
the context of graphene [7J [5], Beenakker pointed out 
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that in addition to the familiar Andreev retro-reflection 
where the reflected hole has a group velocity opposite to 
the incident electron when E < n, there is also the case 
of specular Andreev reflection where the reflected hole's 
group velocity is in the specular direction for E > /i. 
Typical scattering trajectories in these two regimes are 
contrasted in Fig. [T]d) and[lj:). For fx = as considered 
in Ref. [4], the Majorana fermion excitation with lin- 
ear dispersion is associated with the specular Andreev 
reflections in Fig. [TJj). For large [i, as in the case of as 
grown Bi 2 Se3 crystals, one expects very different behav- 
iors at low energies. For the E < /i case, it can be shown 
analytically that the phase of the retro-reflected hole is 
equal to the incident angle of an incoming electron at 
zero energy, 9 = azcsm(hvp k y /ii). This is unique to TIs 
because the wavefunction of a Dirac electron [or hole], 
(f,±e lS ,0,0)[(0,0, f,±e ie )], is determined by the angle 
9, or k y . The resultant hole incident on the opposite S 
with phase of tt retro-reflccts into an electron. This elec- 
tron has exactly the same phase as it started with, thus 
forming an Andreev bound state. 

The remaining key question is whether there will be 
any states at or near zero energy when fi is finite. We can 
answer the question by solving Eq. ([I]) for an idealized, 
step function profile of A(x), 



A(as) = A[6(-x) 



(2) 



The dark lines in Fig. [Tji) shows the zero energy solution 
in the (/x, k v ) plane, with fixed A and the junction length 
w = IOHvf/A.. In general, there exist multiple zero en- 
ergy bound states at discrete k y values {k y } for finite /it. 
For increasing \x and w, these solutions become increas- 
ingly close-packed. This nontrivial result has important 
implications for experiments. The Majorana quantum 
wire is only ideal in the limit of [i, w — > 0. As \x is tuned 
away from the Dirac point, the single zero energy state at 
k = will be replaced by multiple zero energy solutions 
along the k y axis, and eventually a nearly flat dispersion 
at zero energy. 

To unambiguously establish this claim, we solve the 
differential equation / H(x,k y )ip(x,k y ) = Etp(x,k y ) nu- 
merically for a finite size system, x € [0,L] as shown 
in Fig. [TJa), with open boundary conditions at x = 
0, L [9]. Here the quasiparticle wave function tp = 
(u^, U|_, Vf, v^) T , with the label (x,k y ) omitted. To fully 
describe the proximity effect including the induced super- 
conducting correlations in the TI region and the suppres- 
sion of superconductivity near the TI-S boundary, we de- 
termine the order parameter profile A(x) self-consistently 
through the gap equation 

A(x) = g(x) dk v u„,f(x,k v )v* jir (x,ky). (3) 

<L n <LO D 

Here n labels the eigenstates with energy e„, g is the 
effective attractive interaction, and wd is the Debye fre- 
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FIG. 2: (color online) The local spectral function Af(E, k y ,x) 
(upper panel) and local density of states N(E, x) (lower panel, 
red solid line) at the center of the junction, x — 0.5L. One 
sees "flat" Andreev bound states near zero energy for — kw < 
k y < kp, and correspondingly a pronounced peak at zero 
energy in the LDOS in the lower panel. The lower panel also 
shows different LDOS away from the center, for x from 0.52L 
to 0.58L. 



quency. We assume g is zero in the TI region and con- 
stant inside S. We expand ipix, k y ) and A(ir) in Fourier 
series and convert the differential equation into an alge- 
braic equation jlOl III) . Starting with an initial guess 
of A(x) which features phase difference tt, the iter- 
ative procedure is repeated until desired convergence 
is achieved. Note that the phase difference tt is self- 
maintained throughout and not fixed by hand after every 
iteration. Then, the local spectral function, 



A a (E, k y , x) = ^ S ( E ~ e n)\Una(x, k y ) 



and the local density of states (LDOS), 



(4) 



N(E,x) = / dk y J2S(E - 

)r, (5) 

J n,<T 

can be computed for a I. The calculation is checked 
to reproduce known results, e.g., the linearly dispersing 
Majorana spectrum at fj, = predicted in Ref. [1]. 



The upper panel of Fig. [2] shows the spectral func- 
tion at the center of the junction, Af(E,k y ,x = 0.5L) 
(A± is the same for this value of x), with /x=20meV, 
A = 5.5mcV, w = 0.04L, L = 2576nm, /w F =4.1 AeV, 
and the Fermi momentum kp = fi/(Ki)p). In contrast 
to the E ~ Hvpky mode for /j = 0, we see Andreev 
bound states (ABS) near zero energy within a wide re- 
gion — kp < k y < kp, where the slope hv y = dE/dk y 
approaches zero. The appearance of numerous cross- 
ings at exact zero energy for finite fc^also agrees with 
the model calculation above in Fig. mi). Beyond this 
range, e.g. for k y > kp, the spectrum is reminiscent of 
the particle-hole folded dispersion of the helical metal, 
E ~ ±hvp(k y — kp). 

As an approximate ansatz to describe the almost flat 
dispersion, we introduce the following phenomenological 
model for the ABS for large /j> A, 

E/A = c(k/k F ) N , (6) 

where c is a constant and is a large number. To fix 
N, we demand that the slope of the dispersion at energy 
E ~ A coincides with that of the bare dispersion, i.e., 
dE /dk y \p=A = hvp. This gives an estimate of N, 

N ~ fi/A. (7) 

Note that we are only concerned with the ABS dispersion 
near zero energy and its continuation beyond kp. For 
wider junctions, additional subgap ABS appear at finite 
energies, and they are not described by Eq. Q. Our 
ansatz is inspired by the mathematical theory of Dirac 
points with multiple topological charge N as found in 
multilayered system discussed in Ref. [12] . 

The flat dispersion implies a peak at zero energy in 
the local density of states. The lower panel of Fig. [2] 
shows the LDOS at the center of the junction, at the S- 
TI boundary, and slightly into the superconductor for the 
same junction parameters given above. While the zero 
energy peak becomes less pronounced when away from 
the junction center, it remains clearly visible and persists 
even into the superconductor. Thus, the predicted flat 
ABS has a clear experimental signature in the tunneling 
conductance measurements. 

The existence of two regimes including the flat An- 
dreev bound states near zero energy is a general feature. 
We have carried out systematic, self-consistent simula- 
tions for the general case of an inhomogeneous chemi- 
cal potential, e.g., fi(x) = fipi within the TI region and 
/j,(x) = /is 7^ uti inside the superconductors. For ex- 
ample, we fixed us = lOOmeV, A=10mcV, L = 644nm 
and gradually increased /ipi from zero to fj,s to monitor 
the evolution of the local spectrum with respect to the 
Fermi wave mismatch. We see the linear Majorana mode 
gradually changing into the flat ABS. 

Having established the existence of nearly flat ABS 
around zero energy, now we systematically trace the evo- 
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FIG. 3: (color online) Upper panel: Schematic of the peri- 
odic proximity structure with A(x) = A sin(7ra;/a). The wave 
function |w(a;)| for the zero energy states are peaked at the 
domain wall boundaries, x = ma. Lower panel: Energy spec- 
trum for a = 24/wf/A and /i = 4 A is flat at zero energy, 
which has fine structures upon closer inspection. 

lution from the infinitesmal /i, linear dispersing (Majo- 
rana) regime to the large fi flat ABS regime. Also we 
would like to understand the details of ABS within its 
narrow "band width" . To this end, we will consider a 
simple model which generalizes the tt Josephson junction 
to periodic systems. Namely, in Eq. ([!]), the order pa- 
rameter modulates sinusoidally in the x-direction with 
period 2a as schematically shown in the upper panel of 
Fig.§ 

A(x) = Asin(7T2;/a). (8) 

The sign of the order parameter alternates. Thus the 
structure is effectively a periodic array of the tt junctions 
discussed above in the limit w — > 0. One also recog- 
nizes that A(x) describes a stripe/pair density wave, or 
Larkin-Ovchinnikov superconductor |13j . While such su- 
perconductors are hard to find, one may imagine bring- 
ing them in contact with a TI to realize the model con- 
sider here. Now the Hamiltonian H has discrete transla- 
tional symmetry in the x-direction, H(x) — H(x + 2a). 
We can apply the Bloch-Floquet theorem and intro- 
duce quasi-momentum k x living in the Brillouin zone of 
(— 7r/2a, 7r/2a). For the prescribed A(x), the energy spec- 
trum E(k x ,k y ) can be obtained by diagonalizing % in 
fc-space. Note that the TI (non-superconducting) region 
is shrunk to a point, only the homogenous /i is left as 
tuning parameter. 

The lower panel of Fig. [3] shows the spectrum E{k x 
0. k y ) for a = 24hv F /A, fx = 4A. These flat ABS at zero 
energy do not show significant variation with k x . We 
have checked that the wave function of these zero energy 
states are localized at the domain wall boundaries of the 
order parameter field, i.e., at x = ma (red curve in the 
upper panel of Fig. [3]). For example, the wave func- 
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FIG. 4: Fine structures in the energy spectrum of the peri- 
odic proximity structure with fixed a = 12Kvf / 'A and increas- 
ing (i. The linearly dispersing Majorana spectrum at /i = 
splits and develops curvature to eventually become nearly flat 
within (—hp, hp). The number of zero energy crossings in- 
creases with fi. 



tion of the k y = 0, k x = 0, E « mode can be fit well 
with periodic Gaussians \u(x)\ oc exp(— 1.85(7ra/\/2a) 2 )- 
Since a is large in this case, these results agree well with 
the single junction result before. The dispersion, for ex- 
ample, can be fit well using the ansatz in Eq. ([6|. The 
vanishing band width is, of course, only valid on coarse 
scales. Closer inspection, by blowing up the spectrum 
near zero as illustrated in Fig. |4j reveals the busy life 
of the ABS with N c crossings at zero energy, where N c 
scales linearly with fj,, in agreement with Fig. Re- 
markably, all these fine details are compressed within a 
small energy range. 

Fig. [4] illustrates the evolution of the ABS at low ener- 
gies for the periodic structure as fi is increased from zero. 
For small value of fi = 0.83A, the linear Majorana disper- 
sion splits into two, each developing a curvature, as the 
zero energy crossings move to finite k y values. Further 
increasing fi, these two crossings are stretched further 
outwards, while the dispesion within k y € {—kp, kp) be- 
gin being bent and stretched to form the precursor of the 
flat band. At the same time, addition of new crossings 
introduces more twists. The number of crossing scales 
with iV c ~ /t-t/A. The spaghetti now becomes a rope, and 
looking from afar, it appears as a thin thread. 

Flat bands are more novelties than the norm in 
condensed matter [14]. Recently, several authors 
have demonstrated that surface Andreev bound states 
with flat dispersion arise in certain topological su- 
perconductors, for example Cu2;Bi2Se3 [TS] and non- 



centrosymmetric superconductors |161 117] . In addition, 
flat bands within a superconducting gap could precipi- 
tate secondary interactions between quasiparticles. Their 
existence can be traced back to the nontrivial topology 
associated with the gapped bulk, and thus are topologi- 
cally protected. This mechanism giving rise to flat bands, 
via the bulk-boundary correspondence, differs from what 
is considered here. For example, in Ref. |15j . a robust 
crossing at k = is a crucial point in the argument, and 
the total number of zero energy crossings is guaranteed 
an odd number. In our case, states at k y — are gapped 
for finite size systems (or finite period 2a) . Despite these 
differences, the zero modes share the common trait that 
they are associated with the sign change of the order 
parameter when electrons arc reflected at the surface or 
interface. 

Several groups have successfully fabricated Josephson 
structures on Bi2Se3 of various length using a variety 
of superconducting materials including Al, Al/Ti, W, 
Nb, and Pb etc. [T8H22"] . Gate tunable supercurrent 
has been observed and argued to be due to the TI sur- 
face state [TS]. Superconducting quantum interference 
devices based on such junctions have also been demon- 
strated [2TJ [23] . Thus the flat Andreev bound states at 
zero energy, and the zero bias conductance peak in the 
local density of states, predicted here should be experi- 
mentally accessible. Future work will explore control of 
these slowly dispersing Andreev levels working as qubits 
[24] when confinement in the y direction is also intro- 
duced. Our work also suggests the ac dynamics of the 
S-TI-S junctions will likely to be very complex featuring 
different regimes. The flat ABS at zero energy predicted 
for periodic junction arrays may potentially find techno- 
logical applications. For example, a diverging density of 
states at the midgap may be used to generate microwave 
resonances. 

We would like to thank Noah Bray-Ali, Liang Fu, and 
Takuya Kitagawa for helpful discussions. This work is 
supported by ONR grant No. N00014-09-1-1025A. 



[1] M. Z. Hasan and C. L. Kane, Reviews of Modern Physics 
82, 3045 (2010), URL |http://link.aps.org/doi/107| 
1103/RevModPhys.82.3045| 



[2] X. Qi and S. Zhang, Reviews of Modern Physics 83, 



1057 (2011), URL http://link.aps.org/doi/10.1103/ 
RevModPhys. 83 .1057| 
[3J M. I. Katsnelson, K. S. N ovoselov, and A. K. Geim, Nat 
Phys 2, 620 (2006), URL |http : //dx.doi . org/10 . 1038/| 



nphys384 
[4J L. Fu and C 



[ 



L. Kane, Physical Review Letters 
100, 096407 (2008), URL |http://link.aps.org/doi/| 



10 . 1103/PhysRevLett . 100 . 096407 
[5] J. R. Williams, A. J. Bestwick, P. Gallagher, S. S. 
Hong, Y. Cui, A. S. Bleich, J. G. Analytis, I. R. 
Fisher, and D. Goldhaber-Gordon, Physical Review Let- 



5 



ters 109, 56803 (2 012), URL |http: //link .aps.org/ 
doi/10. 1103/PhysRevL ett . 109 . 056803 

[6] J. Moore, Physics 5 (2012). 

[7] C. W. J. Beenakke r, Phys. Rev. Lett. 97, 67007 
(2006), URL |http: //link, aps.org/doi/10. 1103/| 
Phy sRevLett . 97 . 06700?] 

[8] C. W J. Beenakker, Rev. Mod. Phys. 80, 13 37 (2008), 
URL http : //link . aps . org/doi/10 . 1 103/RevModPhys . 
I 80.13371 

[9] The exact choice of boundary conditions will not affect 
the local physics at the junction, by choosing L to be 
large. The open boundary condition may be realized, as 
an example, by having strong ferromagnets at x — 0,L. 
They will intoduce chiral edge modes spatially well seper- 
ated from the center of the junction and will not affect 
the local spectrum there. 

[10] K. Halterman, O. T. Vails, and M. Alidoust, Physical 
Review B 84 , 064509 (2011), URL |http: //link .aps. 

I org/doi/10 .110 3/PhysRevB . 84 . 064509 

[11] M. Lababidi and E Zhao, Phys Rev. B 83, 
184511 (2011), URL |http: //link .aps.org/doi/loT1 

| 1103/PhysRevB.83.184511| ~ 

[12] T. Heikkila and G. Volovik, JETP Letters 93, 59 (2011), 

ISSN 0021-3640. 
[13] A. I. Larkin and Y. N. Ovchinnikov, Sov. Phys. JETP 

20 (1965). 

[14] T. T. Heikkila, N. B. Kopnin, and G. E. Volovik, JETP 

Letters 94, 252 (2011). 
[15] T. H. Hsieh and L. Fu, Physical Review Letters 
108, 107005 (2012), URL jhttp: //link .aps.org/doi7 
| 10 ■ 1103 /PhysRevLett . 108 . 10700_5j 



[16] P. M. R. Brydon, A. P. Schnyd er, and C. Timm, Phys. 

Rev . B 84, 020501 (2 011) , URL |http : //link . aps . org7| 

doi/10 . 1 103/Phy sRevB . 84 . 020501] 
[17] A. P. Schnyder and S. Ryu, Physical Review B 

84, 060504 (201 1) , URL |http : //link . aps . org/doi/10 . | 
| 1103/PhysRevB.84.060504| ~ 

[18] B. Sacepe, J. B. Oostinga, J. Li, A. Ubaldini, N. J. G. 
Couto, E. Giannini, and A. F. Morpurgo, Nat Commun 2 

(2011) , URL j http : //dx . doi . org/10 . 1038/ncomms 1586 
[19] D. Zhang, j! Wang, A. M. DaSilva, J. S. Lee, H. R. 

Gutierrez, M. H. W. Chan, J. Jain, and N. Samarth, 
Phys. Rev. B 84, 165120 (2011), URL|http: //link, aps .| 
org/doi/10 . 1103/PhysRev B . 84 . 165120| ~ 
[20] M. Veldhorst, M. Snelder, M. Hoek, T. Gang, V. K. 
Guduru, X. L. Wang, U. Zeitler, W. G. van der Wiel, 
A. A. Golubov, H. Hilgenkamp , et al., Nat Mater 11, 4 17 

(2012) , URL [ jhttp : //dx . doi . org/ 10 . 1038/nmat3255| 
[21] F. Qu, F. Yang, J. Shen, Y. Ding, J. Chen, Z. Ji, G. Liu, 

J. Fan, X. Jing, C. Yang, et al., Sci. Rep. 2 (2012), URL 
http : //dx . doi . org/10 . 1038/srep00339 
[22] J. R. Williams, A. J. Bestwick, P. Gallagher, S. S. Hong, 
Y. Cui, A. S. Bleich, J. G. Analytis, I. R. Fisher, and 
D. Goldhaber-Gordon, eprint arXiv: 1202.2323 (2012), 
1202.2323. 

[23] M. Veldhorst, C. G. Molenaar, X. L. Wang, 
H. Hilgenkamp, and A. Brinkman, eprint 
arXiv: 11 12.5858 (2011), 1112.5858. 

[24] A. Zazunov, V. S. Shumeiko, E. N. Bratus', J. Lantz, and 
G. Wendin, Phys. Rev. Lett. 9 0, 087003 (20 03), URL 
|http : //link . aps . org/doi/10 . 1 103/PhysRevL ett . 90 . | 

I 087003l ~ 



